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ABSTRACT 

We present the results of a survey of young intermediate mass stars (age 
< 5 Myr, 1.5 < < 15 Mq) in the W5 massive star forming region. We 
use combined optical, near-infrared and Spitzer Space Telescope photometry and 
optical spectroscopy to define a sample of stars of spectral type A and B and 
examine their infrared excess properties. We find objects with infrared excesses 
characteristic of optically thick disks, i.e. Herbig AeBe stars. These stars are 
rare: <1.5% of the entire spectroscopic sample of A and B stars, and absent 
among stars more massive than 2.4 Mq. 7.5% of the A and B stars possess in- 
frared excesses in a variety of morphologies that suggest their disks are in some 
transitional phase between an initial, optically thick accretion state and later 
evolutionary states. We identify four morphological classes based on the wave- 
length dependence of the observed excess emission above theoretical photospheric 
levels: (a) the optically thick disks; (b) disks with an optically thin excess over 
the wavelength range 2 to 24 //m, similar to that shown by Classical Be stars; (c) 
disks that are optically thin in their inner regions based on their infrared excess 
at 2-8 fim and optically thick in their outer regions based on the magnitude of 
the observed excess emission at 24 /im; (d) disks that exhibit empty inner regions 
(no excess emission at A < 8 fim) and some measurable excess emission at 24 fim. 
A sub-class of disks exhibit no significant excess emission at A < 5.8 yum, have 
excess emission only in the Spitzer 8 /zm band and no detection at 24 /im. We 
discuss these spectral energy distribution (SED) types, suggest physical models 
for disks exhibiting these emission patterns and additional observations to test 
these theories. 

Subject headings: accretion, accretion disks — circumstellar matter — infrared: 
stars — planetary systems: protoplanetary disks — stars: pre-main sequence — 
techniques: spectroscopic 
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Introduction 



Recent observations have shown that planets exist around intermediate mass stars (spec- 
tral types A and B, 1.5 < < 15 M©, Rivera et al. 2005, Johnson et al. 2007) as well 
as around solar type (< Mj, ~ 1 Mq) stars. Extensive near-infrared (IR), mid-infrared (IR) 
through submillimeter and millimeter observations of low mass star disks (stellar spectral 
type later than ~K5) have shown that ^90% of these stars lose their optically thic k inner 
disks by an age of 5-7 Myr (jStrom et al.lll989l : iHaisch et al.ll200ll : [Hillenbrand II2008I ). Early 
studies of the more massive Herbig Ae/Be stars (M,^ta r = 2-8 Mp-,) indicate t hat their disks 
may be dissipated within a shorter timescale ~3 Myr (IHernandez et al.ll2005l ). Such surveys 
however, have been constrained to small samples of intermediate mass stars in several nearby 
star forming regions. In order to better compare the low mass star studies with a large sam- 
ple of intermediate mass stars in a single region we have er nbarked on a cornbined infrared 
and optical survey of the massive star forming region W5 ( IWesterhoutI Il958[ iKoenig et al. 
2008h . 



Combined optical, near-IR (JHKs) and Spitzer spectral energy distribution (SED) 
analyses provide us with a simple diagnostic of the dust in young stellar disks. We have 
conducted a study of spectral types and Ha equivalent widths from optical spectra, combined 
with a simple classification scheme for infrared disk SED morphologies of A and B stars in 
W5. § 12] describes our observations. In § 13.31 we present our analysis of the different types 
of disks we detect in the sample of B and A stars in W5, and our procedure for filtering 
foreground and background contamination. In § H] we discuss our results for infrared disk 
types in the context of models of young stellar disk evolution. 



Data 



2.1. Spitzer Photometry 



Our Spitzer IRA C and MIPS obser vations of W5 are described in detail in iKoenig et al. 

fl2008h . The IRAC JPazio et all \20o4 ] observations (PID 20300) were broken down into 
three rectangular Astronomical Observing Requests (AORs) covering ~1.8° x 1.6°, in order 
to observe at multiple rotation angles and help minimize artifacts aligned along columns 
or rows of the array. Each AOR had a coverage of 1 High Dynamic Range (HDR) frame. 
HDR mode results in a 10.4 s and 0.4 s exposure being taken at each position in each map. 
We used the clustergrinder software tools developed by R. Gutermuth to produce final 
image mosaics from these data in each wavelength band. Clustergrinder incorporates 
all necessary image treatment steps, for example, saturated pixel processing and distortion 
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corrections (see iGutermuth et al.l |2008| . for a more complete description of the processing 
performed). Clustergrinder uses the short 0.4 s exposures only in saturated or near- 
saturated regions, so that the combined map has an effective total integration time of 3x10.4 
s = 31.2 s in most of the overlapping areas. We also incorporated in o ur data processin g 
archival data covering AFGL 4029 (from Spitzer GTO program PID 201. lAllen et al.ll2005l ). 
at the eastern end of W5. 



The MIPS ( iRieke et al.ll2004j ) observations were carried out on 2006 February 23 UT 
under our GO-2 program, PID 20300. Images were taken in scan-map mode using the 
medium scan speed for an average exposure time of 41.9 s pixel"^ once frames were combined. 
The raw data were processed with pipeline version S13.2.0. We produced final mosaics using 
the MIPS instrument team Data Analysis Tool, which calibrates the data and app lies a 
distortion correction to each individual exposure before combining (iGordon et al.ll2005[ ). We 
used only the 24 fim band data for our analysis in this paper, since strong background 
emission dominates at the longer wavelength (70 and 160 /im) bands of MIPS, and lower 
sensitivity reduces the number of detectable objects to a level not useful for the present 
study. 

We carried out point source extraction and aperture photometry of all point sources 
on the final IRAC mosaics with Phot Vis version 1.10bet a3. Phot Vis is an IDL GUI-based 
photometry visualization tool (see Gutermuth et al. 20041) that utilize s DAOPHOT modules 
ported to IDL in the IDL Astronomy User's Library ( lLandsmarull993[ ). We used Phot Vis to 
visually inspect the detected sources in IRAC bands 1, 3 and 4, adding sources not detected 
automatically, but clearly visible in the images with the GUI tool and rejected any structured 
nebulosity or cosmic rays mistaken for stellar sources by the automatic detection algorithm. 
To save time, we did not visually check the band 2 photometry in this manner. Instead, 
we took the cleaned band 1 source list as the starting point for finding objects in the image 
and extracted photometry at each position. Radii of the apertures and inner and outer 
limits of the sky annuli were 2.4", 2.4" and 7.2" respectively. The photometry was calibrated 
using large-aperture in-fiight measurements of standard stars, with an appropriate aperture 
correction in each channel to correct for the smaller apertures used in this study. 

Averaged over the whole W5 field, our source catalog is 90% complete to a magnitude 
of 15.5 at 3.6 fim, 15.5 at 4.5 fim, 14.0 at 5.8 /im and 12.7 at 8 fim. The completeness is less 
in re gions of bright diffu se emission: ~14 at 3.6 and 4.5 /im, ~11 at 5.8 fim and ~9.5 at 8.0 
/xm ( IKoenig et al.l 120081 ). 



We conducted point source extraction and aperture photometry of point sources in the 
24 nm MIPS m osaic using the point-spread-function fitting capability in IRAF DAOPHOT 
( IStetsonlll987l ). We visually inspected the image to pick out point sources not automatically 
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detected due to bright diffuse emission evident throughout the image. We match the four- 
band IRAC source hst to the MIPS catalog using a 2" search radius, selecting the object 
closest to the MIPS point spread function centroid in cases where more than one IRAC 
object is a match. 



2.2. Optical Photometry 



We used KeplerCam f ISzentgyorgyi et al.l 120051 ) on the 1.2 m telescope at FLWO to 
image six fields in W5 in Sloan r and z filters on 2006 January 21. KeplerCam has a 
monolithic 4096^ CCD detector giving a 23.1'x23.1' field of view. We binned the images 2 
X 2, giving a scale of 0.68" per binned pixel. Each field was imaged with a series of 10 s and 
180 s exposures to detect both faint and bright objects. Sky subtraction was carried out by 
constructing sky images from median-combined data frames. We extracted simple aperture 
photometry with IRAF DAOPHOT daof ind and phot tasks. 



We observed a further 6 fields in W5 on 2006 September 25 with MegaCam ( iMcLeod et al. 



2000) on the MMT to fill in gaps left by the KeplerCam pointings. MMT/MegaCam has 36 
chips with 2048x4608 pixels, a pixel scale of 0.08" pixel"^ and a total field of view (FOV) of 
24'. We obtained three 0.4 s and three 50 s dithered exposures in the MegaCam r' filter, and 
three 0.4 s and three 70 s dithered exposures in the MegaCam z' filter in gray conditions, with 
0.8"-2.0" image quality in the r' images and 0.96"-1.8" image quality in the z' images. We re- 
duced the data based on the method described in Matt Ashby's MegaCam Reduction Guide. 
Our reduction relied in part on software written specifically for MMT/Mega Cam data reduc- 



tion by Brian McLeod. We used the Two Micron All Sky Survey (2MASS ISkrutskie et al. 



20061 ) stellar catalog to derive p recise astrometri c solutions for each science exposure. We 
used images of the SSA 22 field (iLilly et al.lll99l[ ) taken on the same night to construct sky 
frames for proper sky and fringe subtraction. We also used dithered images of the SSA 22 
field to derive illumination correction images in z' and r' to divide out the variation in zero 
point across MegaCams FOV. We did a weighted coaddition of the reduced images using the 
IRAF imcombine task. We extracted photometry in r' and z' filters with the Phot Vis tool 
as described above. 
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2.3. Optical Spectra 

2.3.1. Hectospec / 6. 5m 
We obtained optical spectra of candidate intermediate mass stars with the Hectospec 



multi fiber spectrograph mounted on the 6.5 m MMT telescope on Mount Hopkins ( iFabricant et al 



19941) • Hectospec is a multi-object spectrograph with 300 fibers that can be placed within a 
1° diameter circular field. We used the 270 groove mm~^ grating, and obtained spectra in 
the range 3700-9000 A with a resolution of 6.2 A. The pointing is fixed using 2 to 3 guide 
stars. Hectospec requires better than 0'.'3 coordinates for fiber positioning; the fibers sub- 
tend a I'.'S diameter circle. We used guide stars drawn from 2MASS; our Spitzer and optical 
photometry coordinates are also registered to 2MASS. Hectospec allowed us to observe a 
large sample in a series of queue observing runs from October 2006 to December 2008. 

Data reduction was performed by S. Tokarz through the CfA Telescope Data Center, 
using IRAF tasks and other customized reduction scripts. The reduction procedure was the 
standard for Hectospec data, with the addition of our special sky subtraction procedure. In 
order to avoid the difficulties encountered when subtracting a sky spectrum made from an 
average over the highly variable H ll region, we took exposures offset by ~5" after each one 
of the configurations. In this way we obtained a sky spectrum very close to each star through 
the same fiber. Background subtraction was performed in IDL to remove each wavelength 
calibrated sky offset spectrum from its corresponding wavelength calibrated object spectrum. 
Since each object spectrum was made up of three averaged observations and the sky spectra 
were only single exposures, a scale factor was calculated at 5600A and 8400A to completely 
remove the sky lines at these points. These two points were used to linearly interpolate a 
scale factor over the whole wavelength range to correct each sky offset spectrum. 



2.3.2. FAST/ 1.5m 



Additional spectra we re obtained using the FAST slit spec trograph mounted on the 
1.5 m telescope at FLWO ( IFabricant Ill994l : IFabricant et al.lll998l ). using a Loral 512x2688 



CCD. We used the standard FAST COMBO configuration, with a 300 groove mm~^ grating 
and a 3" wide slit. This setup gave spectra in the range 3800-7200 A, with a resolution of 
~6 A, comparable to the Hectospec spectra. The stars were observed during several queued 
runs from 2007 September until 2008 December (observers: P. Berlind & M. Calkins). The 
spectra were reduced at the C fA (S. Tokarz) using software developed specifically for the 
FAST COMBO observations (iTokarz fc Roll! 119971 ). A sample of spectra is displayed in 
Figure [H 



4000 5000 6000 7000 8000 9000 

Wavelength (A) 



Fig. 1. — Sample FAST (shorter wavelength scale) and sky-subtracted Hectospec spectra. 
Relative flux (offset by an arbitrary factor) is plotted against wavelength. Spectral types 
assigned by our automated code are given alongside each object. 



- 7- 



2.4. Spectral Classification 



We classified our objects using the spectral classification software SPTCLAS^ (jHernandez et al. 



2004J ). SPTCLASS works well over tlie FAST and Hectospec wavelengtli ranges for stars witli 
spectral types in tfie range from ~09 to ~L0. The main code includes three spectral clas- 
sification schemes: the first is optimized to classify stars in the mass range of T Tauri stars 
(TTS, type K5 or later), the second is optimized to classify stars in the mass range of solar- 
type stars (late F to early K), and the third is optimized to classify stars in the mass range 
of HAeBe stars (F5 or earlier). The schemes are based on 66 spectral indices sensitive to 
changes in Teff but insensitive to reddening, stellar rotation, luminosity class, and S/N. The 
code generates a postscript file which shows the individual results from the different spectral 
indices, and a closer view around the Ha and Li I (6707A) lines. The best result from the 
three schemes was selected for each spectrum by visual inspection — in some cases (about 
5-10% of stars) where the code failed to accurately classify a spectrurn, a type was assigned 



by visual comparison to standards in the literature (jjacoby et al.lll984j : lAndrillat et al.lll995 
Carquillat eraDll997l ). 



3. Analysis 
3.1. A and B Star Sample Selection 



In this paper we specifically target stars of spectral type A and B, roughly corresponding 
to masses 1.5 < < 15 Mq. We select an extinction-limited sample of A and B stars from 
our spectral catalog. Figure |2] shows the observed color- magnitude diagram for the optical 
photometric catalog in W5 (light gray points), overlaid with the photometry for the complete 
spectral catalog in dark gray. We show a theoretical main sequence in r' an d r' — z' (in red) 
derived fror n the absolute magnitude sc ale given in ISchmidt-Kalerl ( 1l982l ) and the optical 
colors from iKenyon fc Hartmanru (119951) conve rted to the Sloan filter set using the color- 
transform relations given in iJordi et al.l (j2006[ ). The main sequence is shown as it would 
appear at a distance of 2 kpc and extinction Ay=^ and Av=4:. We highlight in blue objects 
from the spectral catalog with spectral type A or B whose photometry ±lcr brings them 
within the trapezoidal region marked in the figure. This population of objects makes up our 
sample: a total of 610 spectroscopically confirmed A and B stars. 



http://www. astro. Isa.umich.edu/^hernandj/SPTclass/sptclass. html 
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Fig. 2. — W5 optical photometric sample (light gray points) overlaid with dark gray points 
showing spectral sample. Blue points mark A and B stars in the spectral catalog that fall 
within the extinction limit. A typical error bar is shown in upper right. 



- 9 - 



3.2. The Age of W5 



The full spectroscopic sample in W5 comprises 4800 stars, ranging in spectral type from 
09 to M6.5. Since we don't have a good technique to establish membership o f the r egion, we 
select objects which match the Spitzer excess source list from iKoenig et al.l (|2008[ ): a total 
of 389 stars. We calculate for each star the value of optical extinction Ay from the observed 
optical photometry: 



Av = 2.56 (r' 



[1.584(i?c - Ic) - 0.386]) 



(1) 



We obtain Rq — Ic from lKenyon fc HartmannI (119951 ) according to s pectral type. Equa- 



tion [T] relies on the conversion between Rcic 

and S DSg§ r' z' given in iJordi et al.l (j2006[ ) 
and the extinction versus wavelength values given in Schlegel et al. ( 19981 ) such that Ar' = 
0.843 Av and A,> = 0.453 Ay. 

In Figure [3] (left panel) we present the dereddened optical color-r nagnitude diagram 
for this sub-sample. We overplot isochrones at 1, 5 and 100 Myr from ISiess et al.l (120001 ) . 
assuming a dist ance of 2 kp c to W 5 and using the conversion relations between RcIc and 
SDSS r' z' from Ijordi et al.l (120061 ). The distribution of the later spectral type stars at 
(r' — z')q > 0.3 suggests that W5 is at age 5 Myr or younger. This compares well with 
the age upper limit implied by the presence of the central O stars in the region ( 5 Myr , 
Karr fc Martini l2003l ) and with the bubble expansion ages presented in IVallee et al.l (119791 ) 
(age < 1.4 Myr). However, a trend of increasing stellar age with increasing mass is apparent 
in Fig. 131 In the right panel of the figure we show a zoomed-in portion of the diagram 
focusing on the massive stars, with A and B stars highlighted in blue. In this part of the 
color-magnitude diagram the stars lie fai i iter th an where we would expect them at an age 
of 5 Myr given the models of ISiess et al.l (120001 ). Some small fraction of these objects may 
be background stars, but some are likely intrinsically fainter than the model prediction. 
This result is a co r amoii property of stellar ages derived using the HR diagram as noted by 
Hillenbrand et al.l (120081 ). We believe this does not indicate a true difference in the ages of 
the low and high mass stars in W5. 

The spectral samp le was of sufficient resolution to detec t Lithium A6707 absorption, 
an indicator of youth (IRandich et al.l 1200 ll : iHartmann 1 12003| ) in 15 of these objects. We 
highlight these stars with green boxes in Figure El 



We present the photometry and derived spectral types. Ha equivalent widths (EW) and 



^http://www. sdss.org 
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Ay extinction values for the list of 389 IR excess stars with spectral types described above 
in Table [H Column 14 gives the spectral type and its respective error for each star. The 
uncertainty is a combination of error from the fit of each index to the standard main sequence 
and the error in the measurement of each index. We measured the equivalent width of the 
Ha line using the splot task in the IRAF spectral reduction package noao . onedspec. In 
those cases where both emission and absorption components were visible in the spectrum, 
the EW was measured across the full line profile. The EW(Hq;) is given in column 15 of 
Tab. m We list the extinction Ay in column 16. 



3.3. Disk Classification 



We aim to quantify the relative fractions of different disk morphologies among the A 
and B stars in W5. The infrared spectral energy distributions of young stars can show a 
wide variety of morphologies, depending on the properties and evolutionary stage of the disk. 
In order to display this information, for each star we calculate the excess as a function of 
wavelength by deriving: (J — [X])dered — {J — [X])phot, i-e. the dereddened color at wavelength 
A minus the photospheric contribution at that wavelength. We assume that there is little 
or no excess emission at J band. The photospheric colors for a given spectral type are 
ca lculated by convolving a simple Planck function at the stellar effective temperature given 



m 



Schmidt-Kalerl (119821 ) with the appropriate filter-plus-telescope-plus-instrument response 
curves (as archived for 2MAsdfl and S'jozfeer [2I) . The co lors derived by our bl ackbody 
calculation agree with those in the literature, for example: iBessell fc BrettI (1198811. wit hin 
0.05 mag. We use the infrared extinction law A[\] /A y given by llndebetouw et al.l (120051 ) at 
near-infrared wavelengths and iFlaherty et al.l (120071 ) in the Spitzer bands to deredden our 
infrared photometry for this calculation. We require that objects possess a 3cr excess in at 
least one filter to qualify as an excess object. Upper limit detections were not acceptable as 
evidence of an excess and we also reject objects that are visibly non-point-like in either the 
8 or 24 fim images. 

We identify four morphological classes based on the wavelength dependence of the ob- 
served excess emission above theoretical photospheric levels: (a) optically thick disks that 
resemble the known class of Herbig AeBe stars; (b) disks with an optically thin excess over 
the wavelength range 2 to 24 /xm, similar to that shown by Classical Be stars; (c) disks 



http://www.ipac. caltech.edu/2niass/releases/allsky/doc/ 
"'http://ssc.spitzer.caltech.edu/irac/calib/ 
http://ssc.spitzer.caltech.edu/mips/calib/ 



Table 1. W5 IR Excess Spectral Catalog 
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.61(12) 


G7.0±1.0 
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Note. — This Table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. Values in parentheses by 
photometry signify error in last two digits of magnitude value. Right Ascension and Declination coordinates are J2000.0. Please contact author for full table. 

^ID numbering is the same as that from Koenig et al. (2008). 

^Spectral type and measurement error — we assume that all stars are Dwarf luminosity class V. 

I 
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that are optically thin in their inner regions based on their infrared excess at 2-8 /xm and 
optically thick in their outer regions based on the magnitude of the observed excess emission 
at 24 /im; (d) disks that exhibit empty inner regions (no excess emission at A <8 fim) and 
some measurable excess emission at 24 fim. A sub-class of disks exhibit no significant excess 
emission at A <5.8 /im, have excess emission only in the Spitzer 8 fim band and no detection 
at 24 /im. 

There already exists a plethora of classifications for disks around young stars, in partic- 
ular for low mass stars (M^ <1 Mq) and the classes we have found in the A and B stars in 
W5 should be put into this context. In the low mass case, SEDs are classified according to 
their slope in the near to mid infrared into Class I (an optically thick disk and an envelope). 
Class II (an optically thick disk). Cl ass III (no i i ifrare d excess) and 'transitional disks' (no 



infrared excess shortward of 24 /im). iLada et al.l (120061 ) found a further SED morphology of 



optically thin excesses in the IC 348 cluster that they labeled 'anemic' disks. A and B star 
disks in type (a) resemble the low mass Class I and II disks. Our type (b) disks are similar to 
the anemic disks seen in IC 348. The type (d) disks resemble the 'transitional' disks in low 
mass stars. The thin/thick disks in ty pe (c) are not exactly like any of the se definitions of 



disks, but a similar type was noted by iMalfait. Bogaert. fc Waelkend ( 119981 ) in their survey 
of nearby Herbig AeBe stars. 

In Figure H] we show the excess emission (in magnitudes) over the photospheric level 
as a function of wavelength for four A and B stars in our spectral sample showing excesses 
typical of our new classes. In the list of 610 AB stars in the extinction limited sample 
we found 46 (7.5%) with infrared emission consistent with one of these four morphologies: 
Optically thick/Herbig AeBe (a). Optically thin (b). Optically thin inner and thick outer 
(c) and 'Inner Hole' (d). The 9 optically thick disks constitute a fraction of 1.5% of the 
extinction limited sample. The remaining stars not grouped under the four categories above 
appeared photospheric over the viewable wavelength range although many objects lacked a 
detection at 8 or 24 /im, thus we are likely missing some disks. Inner hole disks are detected 
around stars over the spectral type range A8.5 to BO and mass range 1.6 < < 15 Mq. 

We present the properties of the 46 spectrally confirmed A and B stars with IR excess 
in Table |21 IR excess plots for all of these sources are shown in the Appendix, with the 
letter corresponding to the excess class in upper left of each plot. Disks in the sub-class with 
excess only at 8 /im and no 24 /im detection are marked as 'de' in the plots. 

The breakdown of infrared excess types in our spectral sample is presented in Table [31 
Following the findings of Wolff et al. (2010, in press) we also note a further sub-category 
of group (d) that is: objects with no discernible excess in wavebands up to 5.8 /im and a 
real 3a excess at the 8 /im band. We find three such objects in the W5 sample and mark 



Table 2. AB Stars With IR Excess 





RA 


Dec 


J 


H 




[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


r' 


z 


Sp. Type 


EW(Ha 


) Av 






(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Type'' 


(A) 


(mag) 


Disk': 


9046 


43.495831 


60.666310 


11.00(02) 


10.43(03) 


9.68(02) 


8.49(00) 


8.00(00) 


7.49(00) 


6.93(00) 


5.48(07) 


12.62(12) 


12.24(12) 


B8.5±0.5 


6.42 


2.08(50) 


a 


16953 


45.339983 


60.482399 


11.79(03) 


10.80(03) 


9.87(02) 


8.57(00) 


8.08(00) 


7.44(00) 


6.30(00) 


3.69(07) 


14.04(12) 


13.44(12) 


B8.5±0.5 


-17.82 


2.67(50) 


a 


1562 


41.952386 


60.964075 


12.04(02) 


10.94(02) 


9.96(02) 


8.85(00) 


8.55(00) 


8.08(00) 


7.11(00) 


3.59(07) 


14.06(12) 


13.50(12) 


B9.0±0.5 


-46.50 


2.48(50) 


a 


5084 


42.808983 


60.374804 


11.68(03) 


11.12(03) 


10.44(02) 


9.12(00) 


8.60(00) 


8.19(00) 


7.70(00) 


6.27(05) 


13.23(12) 


12.90(12) 


B9.()±().5 


-0.42 


1.91(50) 


a 


6388 


43.009641 


60.604807 


14.03(04) 


13.35(05) 




10.60(00) 


10.16(00) 


9.73(01) 


8.51(01) 


4.96(05) 


16.25(12) 


15.49(12) 


B9.0±0.5 


-8.40 


3.03(50) 


a 


15405 


44.933304 


60.641602 


13.23(02) 


12.48(03) 


11.60(02) 


10.29(00) 


9.70(00) 


8.97(00) 


7.63(00) 


4.90(04) 


15.11(12) 


14.73(12) 


A2.6±0.6 


-14.19 


1.80(50) 


a 


5232 


42.827700 


60.346598 


12.50(03) 


12.05(04) 


11.60(02) 


10.65(00) 


10.16(00) 


9.77(01) 


9.09(00) 


7.71(06) 


14.15(12) 


13.73(12) 


A6.0±0.5 


5.83 


1.55(50) 


a 


2736 


42.341200 


60.751156 


13.16(03) 


12.11(04) 


11.08(02) 


9.95(00) 


9.55(00) 


9.05(01) 


8.11(02) 


4.23(08) 


15.97(12) 


14.75(12) 


A6.5±0.5 


-11.87 


3.51(50) 


a 


7717 


43.254185 


60.595434 


12.66(02) 


11.99(03) 


11.34(02) 


10.24(00) 


9.81(00) 


9.20(00) 


7.74(00) 


4.71(08) 


14.55(12) 


14.11(12) 


A6.5±0.5 


-2.40 


1.52(50) 


a 


13385 


44.461729 


60.330359 


10.71(04) 


10.46(05) 


10.37(04) 


9.89(00) 


9.71(00) 


9.49(01) 


9.25(02) 


7.25(18) 


11.90(12) 


11.66(12) 


B0.0±0.5 


-12.20 


2.23(50) 


b 


Note. 


— This Table i 


s published 


in its entirety 


in the electronic edition 


of the AstTophysical Journal. A portion is shown here for guidance regarding 


its form and content. Values in parentheses by photometry 


signify error in last two 


digits of magnitude value. 


Right Ascension and Declination coordinates are 


J2000.0. Please contact author for full table. 












^ID numbering is the 


same as that from Koenig et al. (2008) 


. ID numbers greater than 20000 are 


new to this 


paper. 
















^Spectra] type and measiirement 


error — we assume that all stars are Dwarf luminosi 


ty class V. 




















^IR excess morpiiological type as 


described in the text, a: Optically thick; b: optically thin; c: optically thin 


at short A, 


optically thick at long A; d 


: no excess 


at short A, some excess 


at long A. 





I 

I— ' 
CO 

I 



- 14 - 



them as 'de' in Table |2j These may be examples of a disk with a central hole, or objects 
with strong PAH line emission (which would give rise to a sharp increase at 8 fim but not 
at shorter wavelengths). 

We next expand our sample of intermediate mass stars in W5 by adding candid ate A 
and B stars from the full J through 24 fim photometric catalog of iKoenig et al.l (120081 ) that 
were not included in the spectroscopic survey. We choose objects with J-band and IRAC 
channel 1 to 4 photometric uncertainty <0.2, a total of 14283 stars, excluding the list for 
which we have already obtained spectra. As described by Koenig et al. we estimate the 
extinction for each source in the absence of a spectral type from an Ay map constructed 
with 2MASS HKs photometry of sources across W5. We deredden the photometry with the 
value of Ay given by the location of each source within the map. We assume that A and B 
stars at the distance of W5 will have a dereddened J magnitude <13 and use this limit as a 
first cut to find them in the photometric sample: a total of 4805 candidate A and B stars. In 
Figure [5] we show a dereddened plot of 2MASS J — H versus J for the spectroscopic sample 
of A and B stars (with extinction as calculated in § 12.41) and for the photorn etric sample in 
W5 dereddened with the extinction map as described in iKoenig et al.l (120081 ). 



Th e Jn=13 limit rough l y corresponds to a main-seq uence star of spectral type A5 at 
2 kpc (jSchmidt-Kalerl Il982l : iKenyon &: HartmannI Il995l ). With the spectral sample as a 
guide, we examine the infrared colors of sources brighter than this and classify them into the 
same IR SED groups by their location in color space in Figure O We only include objects 
with an excess in at least one color: a total of 64 objects, as before excluding stars for which 
we already have spectra. The full classification scheme is described in the Appendix. Table [3] 
presents the resulting breakdown of disks in the photometric sample by type. The candidate 
photometric A and B stars with IR excess are plotted in Fig. [6] (lower panels). 

The method of finding the photometric-only sample of A and B stars is crude since we 
only require they be bright and possess an apparent infrared excess. As a test, we applied 
the same criteria to our spectroscopic catalog of 4800 objects and found 60 objects. Of these 
stars 62% proved to be A or B type, since some later type stars can be similarly bright and 
possess excess infrared emission. The contamination fraction is not constant between the 
different disk types. The fraction of confirmed A and B stars by group was: Thick: 50%, 
Thin: 63%, Thin/Thick: 45%, Hole: 83%. We add a further row to Table [3] applying these 
success rates to the photometric sample. We note that the proportions of the different disk 
types are consistent within the errors. 
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Fig. 3. — Left panel: dereddened color-magnitude diagram for all Spitzer excess sources 
(Koenig et al. 2008) detected in our spectroscopic survey. Green boxes indicate those 
objec ts with con v incing Lithium absorption. Solid red lines: is ochrones at 1, 5 and 1 Myr 
fror nlSiess et al.l (120001). S olid blue line: Main Sequence from iKenvon fc HartmannI (119951 ) 
and Schmidt-Kalei ( 1982 ): solid green hne: Main Sequence from Siess et al. (j200ol ). Right 
panel: zoomed-in portion showing the A and B stars highlighted with blue points. A typical 
error bar is shown in upper right of each panel. 



Table 3. Disk Type Breakdown 





HAeBe 


Thin 


Thin/Thick 


Hole 


Spectral 


9 


23 


5 


9 


Photometric 


14 


33 


13 


4 


Total 


23 


56 


18 


13 


(%) 


21±4.8 


51±8.4 


16±4.2 


12±3.5 


Alt. TotaP 


16 


44 


11 


12 


(%) 


19±5.3 


53±9.9 


13±4.3 


15±4.5 



^Alternative totals assuming contamination of photomet- 
ric sample by late type stars of 50%, 37%, 55%, 17% respec- 
tively. Errors in percentages are derived assuming Poisson 
statistics. 
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Wavelength X{\im) 



0.5 1.0 
Wavelength X{\im) 



Fig. 4. — Distribution of excess emission in magnitudes: [Aj^ered — [A]phot as a function of 
wavelength for four example stars (remaining objects are shown in the Appendix). Error 
bars include uncertainty in Ay calculation. The zero level represents the photospheric level 
for that spectral type. In upper left are shown the spectral type and I^ (Hq) in A. Black line: 
Herbig AeBe star median and quartiles fr o m iHillenbrand et al I f 19921) . gr ay line: CBe star 



median and quartiles from Cote fc Waters f 198?! ). 
fll99ll . ll994l ). 



Waters et al. 



and [Dougherty et al.l 
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Fig. 5. — Dereddened 2MASS J and H color-magnitude diagram for spectroscopically con- 
firmed A and B stars (blue points) and the photometric sample in W5 (dark gray points). 
Photometry is dereddened for spectroscopically confirmed A and B stars as in § 12.41 and for 
the photometric points with the extinction map of Koenig et al. (2008). Dashed line shows 
the cut we make in the photometric data to isolate candidate A and B stars at the distance 
of W5. 
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Fig. 6. — Upper panels: Color-color diagrams for the spectroscopic sample of A and B stars 
in W5 possessing at least a 3 sigma excess in one or more bands, color-coded according to 
disk class. Red: type (a), blue: type (b), green: type (c), yellow: type (d). Lower panels: 
Color-color diagrams for the photometric sample in W5: gray points show all photometric 
detections, colored crosses show best candidate A and B stars with the same color scheme 
as for the spectroscopic sample. 
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3.4. Accretion Signatures in the AB star sample 

The spectra in this paper have insufficient resolution to precisely measure Ha equivalent 
widths (EW(Ha) and are subject to large errors introduced by our crude sky subtraction 
technique. Thus we are unable to establish reliable measures of the accretion rates in our 
sample of A and B star disks. However, we can still use the measurements of EW(Ha) from 
our spectra to assess to ffist order which of the stars exhibit some noticeable excess emission 
and are still accreting. 

In Figure [7] we plot the measured EW(Ha) as a function of effective temperature for all 
of the A and B stars in the spectral sample. We overlay on this plot the disk excess objects, 
color-coded according to our earlier color scheme. Negative values of EW(Hq;) indicate 
emission lines. 

The Main Sequence level of Ha absorption is clearly traced from ~ 10 A in the A stars 
below 10000 K to ~5 A in the hottest B stars at 30000 K. The largest negative values of 
EW(Ha) (corresponding to emission) are in the thick disk and thin disk sources. In the case 
of the thick disk sources this emission is likely produced by accretion from the circumstellar 
disk. One of the 'thick' objects (ID 5267, a B5.5 star) was observed at multiple epochs (7- 
Oct-2007 and 30-Sept-2008). Its Ha equivalent width varied from -7.99 to 0.49 A between 
these two observations, suggesting variable accretion activity in this source. The 'thin' 
sources' SEDs resemble tho se of Classical Be stars w hich exhibit strong emission in this line 



produced in a stellar wind ( iPorter fc Riviniusll2003l ). The thin/thick and hole sources also 



exhibit some evidence for Ha in emission which suggests that some of these stars may still 
be accreting material from their disks but clearly at a consistently lower rate than in the 
thick disks. 



3.5. Spatial Variation in SED types 



We present the spatial distribution of A and B star disks in W5 in Figure [8]on the Spitzer 
MIPS 24 fim mosaic. We find no systematic trends of location with disk type, however we 
note that the high extinction on the molecular clouds (marked by the bright diffuse 24 /xm 
emission around the rim of W5) wi ll strongly reduce the number of disk excess objects we find 
in our survey in these regions. In iKoenig et al.l ( l2008l ) we calculated the 90% completeness 
limits for our Spitzer survey as a function of location and wavelength. In the H II region 
cavities we find 90% limits of 16.1 (4.5 /im), 14.6 (5.8 /xm), 13.4 (8 fim) and 11.2 (MIPS, 
24 /xm). In the bright diffuse emission regions, the limits are: 14.2 (4.5 /im), 12.3 (5.8 /im), 
9.5 (8 /im) and 8.5 (MIPS, 24 /im). Adopting theoretical colors as calculated in §3.31 these 
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results mean we are able to detect all A and B star excess types in the H ll region cavities. 
On the cloud however, we are limited by a convolution of the detection limits in each band 
with the range of excesses for each disk type. Table S] presents a summary of the resulting 
completeness limits. As a consequence we miss almost all A-type Thin and Hole disk sources 
against the bright cloud diffuse emission. 



Discussion 



We have found a sample of young A and B stars in W5 that possess a variety of disk 
morphologies and that the disks can be grouped into four distinct types. Although this 
scheme is a crude way to describe what is almost certainly a continuum of disk morphologies, 
we can use these types to in vestigate the properties and orig i ns of t hese disks. Indeed a similar 
scheme was constructed by iMalfait. Bogaert. fc WaelkensI (119981 ) from their observations of 
near and mid-infrared excesses in nearby A and B stars. We have been able to note this 
same range of disk types in a single young region. 

We firstly assume that all stars of intermediate mass begin their life in a protostellar 
phase similar to that of low mass stars. As such, they initially possess an accretion disk and 
an infalling envelope. As tin ie progresses, the envelope disap pears and an op t ically thick 
accretion disk is left behind. [Andrews fc Williams! ( 20071 ) and iMuzeroUe et al.l (2000) have 
shown that in low mass stars, the disk-to-star mass ratio can span two orders of magnitude 
(from 0.1 to 10% of the stellar mass) as can the disk accretion rate onto the star. If the same 
property applies to intermediate mass stars, then the initial conditions alone can produce 
very different end states after a given amount of time in the AB stars' disks in a cluster as 
these disks drain onto their stars and dissipate through spreading. 

Viscous draining an d spreading of a disk alo ne are likely too slow to explain observed disk 
evolutionary timescales (IHartmann et al.lll998l ). Several additional mechanisms exist that 
will alter the evolution of the star-disk system in time and operate on more rapid timescales. 
The disk may photoevaporate due to radiation from the host or external stars. Dust grains in 



Table 4. Disk Type Completeness 



Thick Disks Thin Disks Thin/Thick Disks Hole Disks 



90% Completeness Limit A9.5 
Brightest 50% Only N/A 



AO 
Al 



A9.5 

N/A 



AO 
A4.5 
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the disk can settle to the mid-plane and grow in size. Finally, large solid bodies (planetesimals 
and planets) may form in the disk. There is evidence that disk e volution proceeds faster 
in more massive stars than in low mass stars ( iHaisch et al.ll200ll ). [Hernandez et al. J2005h 
suggest that the timescale for removal of the optically thick inner disk in Herbig Ae/Be stars 
(types B5 to A9) is short, given that the fraction of stars with an optically thick inner disk 
as seen in the JHKs bands is lower at a given age amongst A and B stars than in low mass 
GKM stars. In W5 we find no optically thick disks in our spectral sample in stars earlier 
than B8.5. Put another way, only 20±5% of stars with any kind of excess are still optically 
thick in the spectroscopic-only sample, and none above a mass 2.4 Mq. 

External erosion of low r nass star disks has been clearly observed in the Orion Neb- 
ula near the Trapezium stars (lO'Dell et al.lll993l ). but the expected photoevaporation rate 
produced drops off exponentially with distance if the star and disk are beyond 1 pc from 
the powering O star. Photoevaporation of the disk by its host star is driven by heating of 
the disk surface by X-ray (hz/ >0.1 keV), Far-ultraviolet (FUV, hz/ <13.6 eV) and Extreme 
ultraviolet (EUV, hz/ >13.6 eV) radiation. FUV radiation heats the disk. Beyond a certain 
radius in the disk, the gas thermal velocity exceeds the escape velocity of the system and 
material can leave in a photoevaporative wind. This effect removes the outer disk while the 
inner disk drains onto the star through accretion. X-rays and EUV radiation may be able 



to cre ate a gap in the disk and subsequently rapidly clear the entire inner disk ( IClarke et al. 



20011 ). Key to the operation of this mechanism is the ability of the EUV and X-ray photons 



to pe netrate the accretion flow and ionize the disk — a still open question (lErcolano et al. 
20081 ). A disk cleared out from the inside in this manner might resemble the 'inner hole' 



sources in our survey. The t imescale for this process is short however, typically less than 
10^ yr jAlexander et al.lbooel l. 



As a disk evolves, dust grains will settle to the disk mid-plane and can grow in size. 
Grain growth and planetesimal formation will reduce the optical depth in small grains which 
could explain the reduced emission in the inner disk as seen in the thin/thick sources. Gas 
will still be present in such a disk and so the accretion rate should remain at a similar level 
as in the optically thick disks. The formation of a giant planet in the disk (Mf^Mjup) could 
also produce the thin/thick morphology. The planet both reduces the dust optical depth and 
decreases the accretion rate of ma terial onto the st ar by slowing the rate of flow of material 
from the outer to the inner disk ( iRice et al.l 120061 ). In this case the accretion rate will be 
around 10% of the typical thick disk rate. A large enough planet (M3>Mjup) can completely 
clear the inner disk, resulting in a disk more like the 'hole' sources. A study of the disks 
around low mass stars in Taurus appears to show that stars with IR morphologies similar 
to our thin/thick sou rces do indeed have accretion rate s about 10 x lower than the optically 
thick accretion disks ( iNajita. Strom fc MuzeroUd 120071 ). If the same is true for the AB stars 
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in W5, then these sources may be good candidates for giant-planet forming disks. The Ha 
hne in the thin/thick sources are systematically smaller than in the thick sources in Fig. [71 
however these results are not precise enough to establish their relative accretion rates. The 
hole sources are similar to the thin/thick disks in that they show some process that has 
cleared out the inner disk. These could be examples of giant planet formation as well, but 
could also simply be photoevaporated disks. Again their Ha lines show less emission than 
the thick disk sources which could indicate lower accretion rates in these objects. 

The thin disk objects in our sample resemble so-called homologously depleted or weak 
disks. These disks may represent the action of gr ain growth or sonae pro cess that drives mass 



l oss fr om the disk at a wide range of disk radii. iLuhman et al.l (|2010[ ) and iMuzeroUe et al. 



(I2OIOI ) discovered similar disks in their recent studies of low mass young stars, but in far 
fewer numbers. In Figure Owe make a plot similar to their Figures 18-25. They typically 
find <10 objects below the boundary marking the location of Taurus primordial disks (with 
the exception of the sample from IC 348). As a fraction of all disks, our sample of thin disks 
in W5 is larger than in any of their low mass disk samples. We do not know at present why 
they are more prevalent among young intermediate mass stars. The earliest types showing 
this excess are a BO, a BO. 5 and a B2 star. 



5. Conclusions and Future Work 

We have carried out a combined infrared and optical photometric and spectroscopic 
survey of young stars in W5. We present evidence from the larger spectroscopic sample that 
W5 is a young region (age < 5 Myr), in agreement with the presence of massive O stars and 
a short expansion timescale for the Hll region bubbles. More detailed analysis of the HR 
diagram would help to place better constraints on the age of W5. 

Within the spectroscopic sample we identify a population of A and B stars. Of these, 
we find that 46 objects show excess infrared emission in at least one Spitzer band. We find 
that many of these infrared excess sources suggest evolved disks with little or no excess 
infrared emission above the photosphere up to a wavelength 5.8-8 /im, but strong excess at 
longer wavelengths. These may be transitional disk objects that have cleared out their inner 
disks. Photoevaporation is one possible mechanism to create this SED morphology, however 
there is a distinct possibility that their SEDs are produced by the formation of one or more 
planets. If we can assess the contribution of the various mechanisms of disk evolution and 
their frequency of occurence, we can begin to understand what types of planets will form 
in a given disk around a star of given mass, how frequently planet formation occurs in the 
Galaxy and at what point in a star's evolution it begins. A clearer picture of the planet 
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forming potential of intermediate mass stars should soon emerge. 
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Facilities: MMT(Hectospec, Megacam), FLW0:1.5m(FAST), FLW0:1.2m(Keplercam), 
Spitzer (IRAC, MIPS) 



We classify disks from the photometric sample of stars possessing dereddcned J mag- 
nitudes brighter than 13 according to the following scheme using dereddened photometry in 
all cases. Firstly we identify Herbig AeBe star candidates. 



A. 



Disk Classification Scheme 



J - [4.5] > 2.0 
[5.8] - [8.0] > 0.3 



(Al) 
(A2) 



J-H> 




Ks) - 0.6 



(A3) 



J-H> 




Ks) - 0.59 



(A4) 



J-H> 




Ks) - 0.116 



(A5) 



J-H <1.2 



(A6) 



Next we identify the 'thin' disk candidates: 



J - [4.5] > -4.5([5.8] - [8.0]) + 3.0 



(A7) 



- 0.25 < J - [4.5] < 2.0 (A8) 

and either: 

[5.8] - [8.0] < 1.0 (A9) 

or: 

1.2 < [5.8] - [24.0] < 3.0 (AlO) 

We next pick out the 'thin/thick' sources. In some cases this re-classifies thin-disk 
sources. 

0.5 < J - [4.5] < 2.0 (All) 

and either: 

[5.8] - [8.0] > 1.0 (A12) 

or: 



[5.8] - [24.0] > 3.0 (A13) 
Lastly we find the 'hole' sources: 

- 0.25 < J - [4.5] < 0.5 (A14) 

and either: 

[5.8] - [8.0] > 1.0 (A15) 

or: 

[5.8] - [24.0] > 3.0 (A16) 
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B. IR Excess Plots 



In Figure [TO] we display IR excess plots for all of the remaining 42 spectrally confirmed 
A and B stars that belong to groups a, b, c or d. Sources are listed by group, and within 
each of these divisions are sorted by spectral ty pe. Every plot also is ove rlaid with the Group 
I Herbig AeBe star median and quartiles from iHillenbrand et al. (119921 ) , and a line showing 



the Classical Be star r nedian and quartiles from lCote fc WatersI ( 1987 ). Waters et al. ( 1987 ) 
and [Dougherty et al.l ( 1l99ll . ll994j ). where the excesses used to find the medians are calculated 
in the same way as for the W5 stars. 
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Fig. 7. — Ha equivalent width as a function of effective temperature for the A and B stars 
in our sample. Objects showing disk excess are color-coded according to the color scheme 
established in Figure 6. 
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Fig. 8. — SpitzerMlPS 24 /iin mosaic overlaid with AB star disk candidates. Dot points show 
spectroscopically confirmed A and B stars. Cross points mark photometric-only sources. 
Color-coding is as in Fig. El Known O stars in W5 are marked with white asterisks. 
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Fig. 9. — Ks — [S.Q]vs.Ks — [24] color-color plot of spectroscopic A and B stars in W5, 
analogous to Figures 18 to 25 in Luhman et al. (2010). The lower boundary of primordial 
disks in Taurus that they identify is shown as a dotted line. 
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Fig. 10. — IR Excess plots for the spectrally confirmed A and B stars with infrared excess. 
Black points show the data and la error bars. Upper line and e r ror b ars shows Group 
1 Herbig AeBe star median and quartiles from [Hillenbrand et al.l ( 19921). lower l ir ie an d 
error bars sh o ws th e Cl assical Be star ra e dian and qu artiles from ICote fc Waters! (Il987l ). 
Waters et al.l (119871 ) and iDougherty et al.l (Il99ll . Il994j ) . In upper left of each plot is given 



the source ID, spectral type, EW(IIa) and disk type. 
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